PROCEEDINGS 

°F SCIENCE 



0^ 
O 
O 
(N 

a 



Quark mass renormalization with non-exceptional 
momenta 



Yasumichi Aoki* 

RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 



E-mail: yaoki@bnl.gov 

RBC and UKQCD collaborations 



> 

in 

a^ 

(N 

O 

0^ 
O 



Renormalization conditions imposed on quark bilinear vertex functions in the conventional 
RI/MOM scheme use exceptional momentum configurations. With practical values for the lat- 
tice cutoff, these vertex functions are contaminated with unwanted low energy physics (pion pole, 
zero modes, etc), which is a large source of systematic error. These effects can be reduced by us- 
ing non-exceptional momenta. We discuss the quark mass renormalization with non-exceptional 
momenta using 2+1 flavor domain wall fermions based on a recently proposed RI/SMOM scheme. 
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1. Introduction 

Quark masses aie fundamental parameters in QCD. Recent lattice computations made it possi- 
ble to calculate light quark masses up to strange consistently taking into account the light sea quark 
effects. The 2 + 1 flavor domain wall fermion (DWF) calculation reported ||^], 

m^(2GeV) = 3.71(0.16)stat(0.18)syst(0.33)renMeV, (1.1) 
mp(2GeV) = 107.3(4.4)stat(4.9)syst(9.7)renMeV. (1.2) 

The first error is statistical. The second error is the systematic uncertainty in the determination of 
the bare quark mass in the lattice theory. It is dominated by the discretization error and will be 
significantly reduced when the ensemble on a finer lattice is analyzed: we are currently generating 
these configurations which will enable us to perform a continuum extrapolation. In this paper we 
discuss the reduction of the third error which arises in the renormalization of the mass. There 
are two dominant contributions to the error: (i) The non-perturbative renormalization of the mass 
in the RI/MOM scheme. We estimate the corresponding uncertainty to be about 7% due to the 
contamination by chiral symmetry breaking effects as explained below, (ii) The matching from the 
RI/MOM to the MS scheme. The perturbative series for this matching is known to 3-loops ^ 
but converges very poorly. The uncertainty is estimated to be about 6%. 

As was demonstrated in Ref. the unwanted non-perturbative contaminations in RI/MOM 
scheme due to spontaneous chiral symmetry breaking could be reduced by changing the scheme 
to one in which no exceptional momenta are present. The argument stems from the Weinberg's 
theorem [^] on the behavior of the vertex function for large external momenta, where a set of 
external momenta which has zero partial sum is called exceptional. 

In this article after the construction of an RI/MOM scheme with non-exceptional momenta 
for the quark bilinears is summarized, the method is applied to a data set used in the conventional 
RI/MOM renormalization with 2 + 1 flavor DWFs [Q]. The new results are compared with the 
conventional RI/MOM results. 



2. RI/SMOM scheme for quark mass 

2.1 Conventional RI/MOM scheme 

We start briefly summarize the original RI/MOM scheme. A mass renormalization factor is 
completely fixed by introducing the two renormalization conditions on the quark propagator. For 
the conventional RI/MOM scheme, the conditions on the Landau-gauge propagator read 
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(2.2) 



which are imposed at the mass-less point. The renormalized quark propagator and mass are related 
to the bare ones through 



Sr{p)=Z,{^)Sb{p), 



mR = Zm{}i)mB. 



(2.3) 
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Eq. (2.1) determines Z^°'^(/i) which is needed for the Eq. (2.2), which in turn fixes Z^^^{^). In 
the continuum theory, RI/MOM scheme wave function renormaUzation condition Eq. ( |2. l| ) can be 



rewritten in terms of the renormalization condition on the bare amputated Green function n of the 
vector current through the Ward-Takahashi identity as 



^|ki^Tr[7^%(/.)] 



(2.4) 



with Zy = 1. A similar relation applies for the axial vector vertex function, but with a contamination 
of a non-perturbative effect with \/p^ suppression f^. 
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with = 1. The momentum configuration must be stated to fix the renormalization condition. For 



these to give equivalent renormalization condition as Eq. (|2JJ), when the momentum p comes in 
through one fermion line the same p must go out via the other fermion line. This is an exceptional 
momentum configuration {p\ + pi = P — P = 0). It is shown from the Weinberg's theorem that 
the difference of the vector and axial vector vertex amplitude is ~ l/p^ [^. This is consistent with 



the existence of the l/p'^ contamination term in Eq. (2^), which MartinelU et al. derived through 
operator product expansion in Ref. 



On the lattice with DWFs, the use of vector or axial vector vertex function (Eqs. (2.4, 2.5)) 



has an advantage over the quark propagator (Eq. (|2JJ)) in calculating the quark wave function 
renormalization. The derivative with respect to momentum is not practical on the lattice as the 
momenta are quantized on the finite volume lattice. A similar scheme sometimes called as RF, in 
which Eq. ( p] ) is replaced with 
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(2.6) 



is free from the derivative. But, naive implementations of Eq. ( |2.6| ) on the lattice introduce the 
tree revel {pa) error, which is sizable at the momentum range we use [0]. DWFs can utilize the 
conserved axial vector cunent and provide a precise estimate of Z^(= Zy) [|l| of the local currents. 



which in turn allows one to use Eqs. (2^), (2^) with a correction of a Za factor to get Z^*^'^. 



By similar reasons, the mass renormalization should be calculated through bilinear operator 
renormaUzation using the relation Z^ = l/Zj = 1/Zp. In principle, the scalar and pseudoscalar 
renormaUzation factors can be determined at large momenta by imposing the conditions 
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(2.7) 



At finite momenta however, Z5 may differ from Zp due to spontaneous chiral symmetry breaking. 
In particular, for the exceptional momentum case, one needs to subtract the pion pole (for P) or 
double pole (for S with quenching) to remove the divergence of Zp (and Z5 quench) in the chiral 
limit. 
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2.2 RI/SMOM scheme 

Non-exceptional momenta do not have zero partial sum, which suppresses, in the Feynman 
diagram, the small momentum flow leading to non-perturbative contamination. Among various 
choices of non-exceptional momenta, we adopt the symmetric one ji^ = p\= p\ = where q = 
Pi —P2- This choice is convenient because only one invariant is involved {2pi ■ p2 is also jji^). Here 
we briefly review the SMOM (symmetric MOM) scheme mass renormalization which makes use 
of the symmetric momentum configuration. The SMOM scheme is discussed in detail in Ref. [||]. 
We will demonstrate in the next section how the use of this renormalization scheme reduces the 
unwanted non-perturbative contaminations compared to the conventional MOM scheme. 

Other than changing the momentum configuration the SMOM scheme follows the same steps 
as MOM scheme. The renormalization conditions are defined using trace conditions with speci- 
fied projectors on the vertex functions YIq. For the scalar and pseudoscalar operators, the same 
projection operators as MOM scheme ^1, ^75 as shown in Eq. (2^) are used. The vector and ax- 



ialvector operator will be used to calculate Zq. Original MOM scheme uses ^Y^, by which 



one can relate the traced vertex function to z^^^ Eq. (lA). If we used these projection operators 
for the symmetric (non-exceptional) momenta, the resulting Zq would completely differ from that 
of MOM scheme (or RF scheme). Instead, we adopt -[^4^1^^ j^Ys^^i^- show that the 

use of these projection operators give Zq through the vector and axial vector Ward-Takahashi 
identities. The matching of Zq of RF and MS has been calculated to three loops [^, ^, which one 
can just use or can use for the check against the calculation with the vertex functions of vector and 
axialvector in the SMOM scheme. It is worth mentioning that RF and RI/MOM Zq are same up to 
one loop. Thus, the resulting Zq from the vector and axialvector current of SMOM scheme should 
be close to that of original MOM scheme. The perturbative matching of quark mass from SMOM 
to MS, m^^{iJ,) = C„,(sMOM Ms) • has been calculated to one loop ^ as 

C,„(sMOM^Ms) = \--^Cf X (0.484 -0.172(§) + O(a^), (2.8) 

where ^ is the gauge parameter. The same quantity for the original RI/MOM to MS has much 
larger correction (both constant and linear coefficient of c^): 

CJmom ^ ms) = 1 - ^Cf X (4 - + 0(af). (2.9) 

471 

The one loop correction for the Landau gauge = 0) is 1.5% for SMOM and 12% for MOM at 
jj. =2 GeV. The three loop correction is still large: 6% for MOM, which was taken as a conservative 
estimate of systematic error of perturbative matching [Q]. The small correction of SMOM scheme 
is realized through cancellation of finite terms depending on the momentum structure. We have not 
understood if this is an universal property with SMOM scheme, which would persists beyond one 
loop. 

3. Numerical test of the RI/SMOM scheme 

We test the RI/SMOM scheme using the A^j = 2 + 1 DWF data set at a" ^ ~ 1 .7 GeV on 
16^ X 32 lattice with Lg = 16, M5 = 1.8. The quark propagators have been calculated with the point 
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Figure 1: — Ay in the chiral limit as function of for MOM and SMOM scheme, where results with 
linear and quadratic chiral extrapolations in quark mass are shown. 



source . 

Let us first look at the difference of tlie vector and axial vector vertex amplitude, 

^MOM ^ ^Tr[7,nvJ, = ^Tr[757Mn^J O-D 

for the exceptional momentum and for the symmetric (non-exceptional) momentum with 

^SMOM ^ _l_Tr[^g^nvJ, Al™ = ^Tr[75^^^n^J. (3.2) 

Fig. |l| shows the differences in the chiral limit as functions of (= q^). The original MOM 
scheme has non-zero difference due to the spontaneous chiral symmetry breaking. The decrease of 
the difference as momentum increases is due to the recovery of the symmetry and the size is up to 
1% of the average (A^ + Ay) /2 in the region of the momentum we use (pa)^ > 1.3. The difference 
is much suppressed for the SMOM scheme. The linear chiral extrapolation gives results consistent 
with zero. The quadratic extrapolation gives non-zero value, but one order of magnitude smaller 
than MOM. 

The scalar and pseudoscalar vertex amplitudes with MOM and SMOM scheme are shown in 
Fig. A large difference between P and S is observed for the MOM scheme. Since Ap diverges 
as ~ 1 /m in the chiral limit due to the existence of pion pole (^\lf) / m\ and A^ stays finite, the 
difference becomes infinite. One should note that these quantities are identical to all order in 
perturbation theory. The symmetry is badly broken for the non-perturbative renormalization. On 
the other hand, Ap and A5 are consistent with each other for the SMOM case at larger momentum 
{paf > 1. 

As our gauge ensembles have been sampled at single value of nearly physical strange mass, 
we have a systematic error from 7^ even after the two-flavor unitary chiral extrapolation. This 
error in the MOM scheme, which turned out to be 7%, was estimated from the response of the As 
to the u, d quark mass in Ref. [Q]. This error may represent the tolerance of this particular quantity 

'The statistical error could be much improved if the volume source was used |H llQ, nll|. 
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Figure 2: and Ap for MOM and SMOM scheme, ntf refers to the average u, d mass. Points connected 
with the solid lines show the values in the two-flavor unitary chiral limit {mf + m^cs ^ 0). 
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Figure 3: Preliminary results of Z,„ '-' (;?) as a function of renormalization scale p and Z,„ (/x ~ 2GeV) 
as a function of SMOM^ MS matching scale p. The extrapolation {pa)^ is shown with statistical error 
only, which is consistent with the (fi = 2GeV) through MOM scheme with the large systematic error. 



to the emergence of the low energy scale (~ Aqcd)- So even if the — > limit was performed, 
error of similar size would remain due to the non-perturbative effect whose energy scale is about 
the same. Now, in the SMOM case, as shown in the figure, mass dependence is greatly reduced. If 
we adopt the same method, the systematic error of SMOM is about 3% for the scalar or negligible 
(comparable to the statistical error) for the pseudoscalar. 

Black symbols in Fig. |3| show preliminary results of Z,,, = (A5 + Ap) / {Za{Aa + Ay)} with 
SMOM scheme (in the chiral limit) as a function of renormalization scale p, where the Za [p 
estimated form hadronic two point functions is used. Matching to MS with Eq. ( p^ and running 
to /X = 2 GeV with the two-loop anomalous dimension, one obtains the blue symbols as a function 
of matching scale p. The extrapolation (pa)^ —>■ using the points [pa)^ ^1.2 gives a consistent 
result with the same quantity but through original MOM scheme [Q] with 3-loop matching and 
4-loop running. 
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4. Conclusion 

The RI/SMOM scheme, constructed in the framework of the conventional RI/MOM scheme 
with the use of non-exceptional momenta, works very well for reducing non-perturbative contam- 
ination for the quark mass renormalization. The systematic error is reduced to 3% level for Z^, 
while it was 7% for the original MOM scheme. This shows the success of the SMOM scheme 
which was designed to reduce the unwanted non-perturbative contamination. Another systematic 
error is from truncation in the perturbative matching to MS. If we estimate the systematic error for 
the SMOM scheme from the size of 0(o;.s ) at our typical momentum size /i = 2 GeV, it is 1.5 %, 
which is much smaller than 6 % at 0{a^) for the MOM scheme. Further discussions are needed 
for the better understanding of the systematic error of the perturbative matching. 

The first non-trivial test of the SMOM scheme was successful. Application to other bilinear 
operators such as tensors would be straightforward. Similar scheme can be constructed for four- 
quark operators for — K^^ mixing in the standard model and beyond, and for K ^ nn decays. 



We thank Christian Sturm for collaborating on the RI/SMOM renormalization. Presented 
numerical data are obtained through reanalyzing the published data computed on the QCDOC 
machines at RIKEN BNL Research Center, Columbia University and University of Edinburgh. 
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